Tensor Forces and the Ground-State Structure of Nuclei 
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Two-nucleon momentum distributions are calculated for the ground states of nuclei with mass 
number A < 8, using variational Monte Carlo wave functions derived from a realistic Hamiltonian 
with two- and three-nucleon potentials. The momentum distribution of np pairs is found to be much 
larger than that of pp pairs for values of the relative momentum in the range (300-600) MeV/c and 
vanishing total momentum. This order of magnitude difference is seen in all nuclei considered and 
has a universal character originating from the tensor components present in any realistic nucleon- 
nucleon potential. The correlations induced by the tensor force strongly influence the structure of 
np pairs, which are predominantly in deuteron-like states, while they are ineffective for pp pairs, 
which are mostly in^So states. These features should be easily observable in two-nucleon knock-out 
processes, such as A{e,e'np) and A{e,e'pp). 

PACS numbers: 21.60.-n,21.30.Fe,25.30.-c 



The two preeminent features of the nucleon-nucleon 
(NN) interaction are its short-range repulsion and 
intermediate- to long-range tensor character. These in- 
duce strong spatial-spin-isospin NN correlations, which 
leave their imprint on the structure of ground- and 
excited-state wave functions. Several nuclear properties 
reflect the presence of these features. For example, the 
two-nucleon density distributions p^g (r) in states with 
pair spin 5'=! and isospin T=0 are very small at small 
inter-nucleon separation r and exhibit strong anisotropics 
depending on the spin projection Ms [H. Nucleon mo- 
mentum distributions N{k) 0, Q and spectral functions 
S{k, E) have large high-momentum and, in the case of 
S{k, E), high-energy components, which are produced by 
short-range and tensor correlations. The latter also in- 
fluence the distribution of strength in response functions 
R{k,uj), which characterize the response of the nucleus 
to a spin-isospin disturbance injecting momentum k and 



energy lo into the system js], Q. Lastly, calculations of 
low-energy spectra in light nuclei (up to mass number 
A=1Q) have demonstrated that tensor forces play a cru- 
cial role in reproducing the observed ordering of the levels 
and, inparticular, the observed absence of stable A = % 
nuclei Q. 

In the present study we show that tensor correlations 
also impact strongly the momentum distributions of NN 
pairs in the ground state of a nucleus and, in particular, 
that they lead to large differences in the np versus pp dis- 
tributions at moderate values of the relative momentum 
in the pair. These differences should be observable in 
two-nucleon knock-out processes, such as A{e, e'np) and 
A{e, e'pp) reactions. 

The probability of finding two nucleons with relative 
momentum q and total momentum Q in isospin state 
TMt in the ground state of a nucleus is proportional to 
the density 



PTMT(q,Q) = 2{2J+l) ^ J dridr2dr3- ■ -drAdr'j^dr'^ Vjmj ^a, ra, ■ ■ • , r^) 



(1) 



where ri2 = ri — r2, R12 = (ri -I- r2)/2, and similarly for 
r'12 and Ri2- ^Vmt (12) is the isospin projection opera- 
tor, and ijjjMj denotes the nuclear wave function in spin 
and spin-projection state JMj. The normalization is 

/ (S^s PT'^i^ = ^tMt , (2) 

where Ntmt is the number of AW pairs in state TMt- 
Obviously, integrating ptmt (Qi Q) over only Q gives the 



probability of finding two nucleons with relative momen- 
tum q, regardless of their pair momentum Q (and vice- 
versa) . 

The present study of two-nucleon momentum distribu- 
tions in light nuclei (up to A=8) is based on variational 
Monte Carlo (VMC) wave functions, derived from a re- 
alistic Hamiltonian consisting of the Argonne wis two- 
nucleon and Urbana-IX three-nucleon [§] interactions 
(AV18/UIX). The high accuracy of the VMC wave func- 
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tions is well documented (see Refs. and references 

therein), as is the quality of the AV18/UIX Hamiltonian 
in quantitatively accounting for a wide variety of light 
nuclei properties, such as elastic and inelastic electro- 
magnetic form factors jl2i] . and low-energy capture re- 
actions However, it is important to stress that the 
large effect of tensor correlations on two-nucleon momen- 
tum distributions and the resulting isospin dependence of 
the latter remain valid, even if one uses a semi-realistic 
Hamiltonian model. This will be shown explicitly below. 



The double Fourier transform in Eq. ([T]) is computed 
by Monte Carlo (MC) integration. A standard Metropo- 
lis walk, guided by | -0 jm, (ri, r2, ra, . . . ,rA)|^, is used to 
sample configurations . For each configuration a two- 
dimensional grid of Gauss-Legendre points, Xi and Xj, is 
used to compute the Fourier transform. Instead of just 
moving the ip' position (r'12 and R'12) away from a fixed 
'0 position (ri2 and R12), both positions are moved sym- 
metrically away from ri2 and R12, so Eq. ([T]) becomes 



PTMT(q, Q) = 2(2J + ^l) ^ J '^^^ ^^'^ dr^--- dvA dxdX 0^^,^^ (ri2+x/2, R12 + X/2, rg, . . . , r^) 



X e 



— iq-x — iQ-X 



Ptm^ (12) ^jMj (ri2 -x/2, R12 -X/2, rg, . . . , r^) 



(3) 



Here the polar angles of the x and X grids are also 
sampled by MC integration, with one sample per pair. 
This procedure is similar to that adopted most recently 
in studies of the ^Ke{e,e'p)d and '*He(e, e'p)'^H reac- 
tions 14], and has the advantage of very substantially re- 
ducing the statistical errors originating from the rapidly 
oscillating nature of the integrand for large values of q 
and Q. Indeed, earlier calculations of nucleon and cluster 
momentum distributions in few-nucleon systems, which 
were carried out by direct MC integration over all coordi- 
nates, were very noisy for momenta beyond 2 fm~^, even 
when the random walk consisted of a very large number 
of configurations 

The present method is, however, computationally in- 
tensive, because complete Gaussian integrations have to 
be performed for each of the configurations sampled in 
the random walk. The large range of values of x and X 
required to obtain converged results, especially for ^He, 
require fairly large numbers of points; we used grids of 
up to 96 and 80 points for x and X, respectively. We 
also sum over all pairs instead of just pair 12. 

The np and pp momentum distributions in ^He, ^He, 
^Li, and *Be nuclei are shown in Fig.[T]as functions of the 
relative momentum q at fixed total pair momentum Q=0, 
corresponding to nucleons moving back to back. The 
statistical errors due to the Monte Carlo integration are 
displayed only for the pp pairs; they are negligibly small 
for the np pairs. The striking features seen in all cases 
are: i) the momentum distribution of np pairs is much 
larger than that of pp pairs for relative momenta in the 
range 1.5-3.0 fm~^, and ii) for the helium and lithium 
isotopes the node in the pp momentum distribution is 
absent in the np one, which instead exhibits a change of 
slope at a characteristic value of p ~ 1.5 fm~^. The nodal 
structure is much less prominent in ^Be. At small val- 
ues of q the ratios of np to pp momentum distributions 
are closer to those of np to pp pair numbers, which in 



■^He, ''He, ^Li, and ®Be are respectively 2, 4, 3, and 8/3. 
Note that the np momentum distribution is given by the 
linear combination ptMt=io+PtMt=oo^ while the pp mo- 
mentum distribution corresponds to ptMt=ii- The wave 
functions utilized in the present study are eigenstates of 
total isospin (1/2 for ^He, and for "^He, '^Li, and ^Be), 
so the small effects of isospin-symmetry-breaking inter- 
actions are ignored. As a result, in ''He, ^Li, and *Be 
the ptMt is independent of the isospin projection and, 
in particular, the pp and T = \ np momentum distribu- 
tions are the same. 

The excess strength in the np momentum distribution 
is due to the strong correlations induced by tensor com- 
ponents in the underlying NN potential. For Q=0, the 
pair and residual {A-2) system are in a relative S-wave. 
In ^He and ^He with uncorrelated wave functions, 3/4 of 
the np pairs are in deuteron-like T, S'=0,1 states, while 
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FIG. 1: (Color online) The np (lines) and pp (symbols) mo- 
mentum distributions in various nuclei as functions of the 
relative momentum q at vanishing total pair momentum Q. 
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the pp, nn and remaining 1/4 of np pairs are in T, 5=1,0 
(quasi-bound) states. When multi-body tensor correla- 
tions are taken into account, 10-15% of the T, 5*^1,0 
pairs are spin-flipped to T, 5=1,1 pairs, but the num- 
ber of T, S'=0,1 pairs hardly changes [l|. In A > A nu- 
clei, some np and pp pairs will be in relative P-waves 
(T, 5=1,1 and 0,0) when one particle is in the s-shell and 
one in the p-shell. Nevertheless, 5.5 out of 9 np pairs in 
^Li (9 out of 16 in ^Be) are expected to be deuteron-like, 
while half the pp pairs will be in T, 5=1,0 states and half 
in T, 5=1,1 states [15i]. The tensor force vanishes in the 
former and is weak in the latter. 

These expectations are born out by our calculations, 
as Fig. [T] clearly demonstrates. The np momentum dis- 
tributions for q values larger than 1.5 fm^^ only differ 
by a scaling factor, and indeed all scale relative to the 
deuteron momentum distribution as shown in Fig. [51 The 
deuteron Pnp{q, Q = 0) has been scaled to match that 
of ''He at q = 1.5 fm"'^. (The scaling property does 
not extend to the low q region, because there binding ef- 
fects take over.) The S- and D-wave components of the 
deuteron density are also shown in the figure as dotted 
lines; the D-state component is the dominant part over 
the range 1.4-4.0 fm~^, while the S-state component has 
a node at 2.1 fm~^. 

Similar considerations in relation to scaling and bind- 
ing effects also remain valid when considering the pp mo- 
mentum distributions. In particular, the node seen in 
helium and lithium is reminiscent of the node in the S- 
wave momentum distribution (shown in Fig. [2]). In larger 
nuclei the node is filled in, due to the fact that pp pairs 
are not exclusively in quasi-bound states, but can also be 
in P-wave or higher partial-wave states. 

Figure [2] also shows the np and pp momentum distri- 
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FIG. 2: (Color online) The np (lines) and pp (symbols) mo- 
mentum distributions in *He obtained with different Hamil- 
tonians. Also shown is the scaled momentum distribution for 
the AV18 deuteron; its separate S- and D-wave components 
are shown by dotted lines. 



butions in ^He obtained with Hamiltonians of decreasing 
sophistication, ranging from the fully realistic AV18/UIX 
model to the semi-realistic Argonne v'q (AV6') to the rel- 
atively simple Argonne U4 (AV4'). The AV4' and AV6' 
potentials are constructed to preserve as many features 
of NN scattering and deuteron properties as feasible 0] . 
The AV4' has only central, spin, and isospin operators 
with no tensor component: it reproduces the ^So phase 
shift, and the deuteron binding energy but with only an 
S-state component; the D-state, induced by the tensor 
term in the potential and associated with one-pion ex- 
change at long range, is absent. On the other hand, 
the AV6' includes tensor terms and produces a bound 
deuteron quite close to that of AV18. These features of 
the underlying NN potential, or the lack of them, are re- 
flected in the calculated momentum distributions. In par- 
ticular, note the node which develops in the AV4' np mo- 
mentum distribution, due to the purely S-wave nature of 
the deuteron-like state. On the other hand, the AV6' and 
AV18/UIX results are very close to each other, demon- 
strating the essential role played by the tensor potential 
in substantially increasing the intermediate-momentum 
components of np pairs. 
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FIG. 3: (Color online) The momentum distributions pNN{q) 
and pnn{Q) in '*He for np and pp. The inset shows the ratios 
Pnp{q)/ppp{q) and pupiQ) / Ppp{Q) ■ 

The momentum distributions pnn{q) and pnn{Q) ob- 
tained by integrating over Q or q respectively are plot- 
ted in Fig. [3] for ''He with the AV18/UIX Hamiltonian. 
There is no node in Ppp{q)'. it is filled in by the con- 
tributions of pp channels other than '^Sq. These chan- 
nels are now allowed, because the orbital angular mo- 
mentum between the pp pair and the residual nn cluster 
is not constrained to vanish for Q > Q. A remnant of 
the strong tensor correlations affecting PnpiQ) still per- 
sists. This is illustrated in the inset of Fig. [31 which 
shows the ratio Pnpi<l) / Pppil) ■ However, the effect is far 
less dramatic than in the back-to-back ((5=0) kinematics. 



4 



Note that Pnn{q) and pnn{Q) have the same normahza- 
tion, i.e. {A - Z)Z for np and Z{Z - l)/2 for pp. The 
PnpiQ) / PppiQ) ratio is close to the np to pp pair number 
ratio — 4 in "^He — over the whole range of Q values. This 
holds true for all the nuclei studied. 

The most direct evidence for tensor correlations in nu- 
clei comes from measurements of the deuteron structure 
functions and tensor polarization by elastic electron scat- 
tering (l6| . In essence, these measurements have mapped 
out the Fourier transforms of the charge densities of the 
deuteron in states with spin projections ±1 and 0, show- 
ing that they are very different. In other processes, such 
as ^H((i, 7)'*He [l^ at very low energy, or proton knock- 
out from a polarized deuteron [13] (as well as in the nu- 
clear properties mentioned at the beginning of this let- 
ter), the effects of tensor correlations are more subtle and 
their presence is not easily isolated in the experimental 
data. This is because of "contaminations" from initial 
or final state interactions and many-body terms in the 
transition operators. 

Some of these corrections will also pollute the cross 
sections for (e,e'np) and (e,e'pp), or {p,pp) and {p,ppn), 
knock-out processes in back-to-back kinematics. How- 
ever, one would expect the contributions due to final 
state interactions in the np and pp reactions, both be- 
tween the nucleons in the pair and between these and the 
nucleons in the residual {A — 2) system, to be of similar 
magnitude for relative momenta in the range (300-600) 
MeV/c. In the electrodisintegration processes, the lead- 
ing electromagnetic two-body currents associated with 
pion and p-meson exchange, denoted respectively as PS 
and V in Ref. |19l] , vanish in pp because of their isospin 
structure. Of course, they will contribute in np, but are 
not expected to produce large effects. Thus the ratio of 
np to pp cross sections should be much larger than unity 
for relative momenta within (300-600) MeV/c, refiecting 
the large difference between the corresponding momen- 
tum distributions in this range. There are strong indica- 
tions from a recent analysis of a BNL experiment, which 
measured cross sections for (p, pp) and {p, ppn) processes 
on ^^C in kinematics close to two nucleons being ejected 
back to back, that this is indeed the case [23| • The ratio 
of {p,ppn) to {p,pp) events over the range of relative mo- 
menta (275-550) MeV/c is found to be roughly 20, albeit 
with a rather large error. Hopefully, a more precise value 
for this ratio will become available in the near future, 
when the analysis of ^^C{e, e'np) and ^^C{e,e'pp) data, 
taken at Jefferson Lab, is completed (2]| . 

It would be interesting to extend these measurements 
to other nuclei. In ■^He and ''He, one would expect the 
node in the pp momentum distribution to be filled in 
by interaction effects in the final state ,14]. However, 
the ratio of np to pp cross sections in the range (300- 
600) MeV/c should still reflect the large value of the np 



momentum distribution at these values of relative mo- 
menta. This would provide a further, direct verification 
of the crucial role that the tensor force plays in shaping 
the short-range structure of nuclei. 

We dedicate this paper to the memory of Vijay R. 
Pandharipande, a mentor and friend, who was deeply 
interested in evidence for correlations in nuclei. This is 
also in remembrance of our colleague and friend Adelchi 
Fabrocini, who contributed significantly to the theoreti- 
cal study of correlations in strongly interacting systems. 
An interesting and stimulating conversation with D. Hig- 
inbotham and E. Piasetzky is gratefully acknowledged by 
one of the authors (R.S.). This work is supported by the 
U.S. Department of Energy, Office of Nuclear Physics, 
under contracts DE-AC05-06OR23177 (R.S.), DE-AC- 
02-06CH11357 (S.C.P. and R.B.W.), and W-7405-ENG- 
36 (J.C.). The calculations were made at Argonne's Lab- 
oratory Computing Resource Center. 
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